On the theory of domain switching kinetics in ferroelectric thin films. 
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We investigate theoretically the polarization switching kinetics in ferroelectric thin films. In 
such substances, the domain walls are pinned by (usually dipole) defects, which are present also 
in ordered samples as technologically unavoidable impurities. This random interaction with dipole 
pinning centers results, in particular, in exponentially broad distribution of switching times. Under 
supposition of low pinning centers concentration, we derive the distribution function of switching 
times showing that it is not simply Lorentzian (as it was first suggested by Tagantsev et al [Phys. 
Rev. B , 66, 214109 (2002)] ), but is a "square of Lorentzian", which is due to the vector nature of 
electric field. This improved formalism delivers a better description of available experimental data 
and elucidates the physical mechanism of domain switching times distribution. 
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Fundamental understanding of polarization switch- 
ing kinetics in ferroelectric thin films is very impor- 
tant as there are numerous applications of this phe- 
nomenon, the most important being in ferroelectric mem- 
ories In bulk ferroelectric crystals, the successful 
description of switching kinetics had been achieved us- 
ing Kolmogorov-Avrami-Ishibashi (KAI) model fo- 
cused on the statistics of nucleation of reversed polariza- 
tion droplets within domains. It has been shown in Rcf. 
Q that application of KAI model to ferroelectric thin 
films does not adequately describe the experiment in the 
case when wide range of switching times is considered. 
Tagantsev et al 0] have shown that this is primarily be- 
cause the KAI model does not account for broad switch- 
ing times distribution, which is especially important in 
the second step of reversed polarization droplets evolu- 
tion, namely the coalescence of initial droplets into larger 
(inversely polarized) domains. This is because in KAI 
model any single droplet of reverse polarization, expand- 
ing unrestrictedly, can provide switching of the whole 
sample. To improve the description of experimental sit- 
uation, Tagantsev et al suggested nucleation-limited 
switching (NLS) model which considers a sample as an 
ensemble of elementary regions which incorporates natu- 
rally the broad distribution of their switching times. 

While in ceramic or polycrystalline samples this dis- 
tribution is due to disorder, in crystalline samples there 
are many unavoidable defects, which also yield the dis- 
tribution of switching times. It is firmly established 
(see [1] and refs therein), that defects and impurities 
in ferroelectrics play a role of pinning centers for do- 
main walls so that under external electric field the do- 
main walls motion has a character of creep. This creep 
motion implies following dependence of a domain wall 
speed on the local electric field £ at a pinning site 
v oc 1/t oc exp[— (U / k,BT)(E / EY], where r is a domain 
wall switching time, U is an energy barrier between dif- 
ferent polarization directions and Eq is a critical field for 
a domain wall creep @, 0] . 



It has been shown (see, e.g. [8[ and refs therein), that 
dynamical exponent fi (defining the details of a domain 
wall pinning potential) almost equals to one so that one 
can get following expression for field dependent domain 
wall switching time 
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where A = V E^jk^T , T m is the switching time at infinite 
local field E. Here, following Ref. we identify the 

local field at a pinning site as a sum of external electric 
field E cx t and internal random field E (below we will see 
that this is actually a modulus of electric field vector) 
which appears due to disorder. 

In non-KAI models of switching, where a distribution 
of switching times is assumed, the expression for the vol- 
ume fraction p(t) of a ferroelectric switched by time t 
assumes the form 0, 
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where y = t/too, y min and ?/ max are, respectively, the 
lower and upper limits of the spectrum of dimensionless 
switching (waiting) times in the system. Their precise 
form will be determined below. As creep character of 
domain walls motion implies broad spectrum of switch- 
ing times, their distribution function <7(log(r/r 00 )) is as- 
sumed to be broad enough to describe this spectrum. Its 
specific form had been chosen as Lorentzian 0, [l(J . 

In this Letter, we show that for low pinning centers 
concentration, when they are spatially uncorrelated (see 
Refs. [ll], [Hf for details), it is possible to calculate 
exactly the distribution function /(E) of local electric 
fields E at the domain pinning sites. We show that this 
distribution function reduces to the one of the electric 
field modulus E = |E| and has the form of "square of 
Lorentzian" rather then simple Lorentzian. This is be- 
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cause pure Lorentzian gives the distribution of a scalar 
quantity, which in the case of above dipole defects ensem- 
ble corresponds to their resonant frequencies flll.[l2j . At 
the same time, the application of formalism [1 21 ] to the 
vector quantity like electric field gives for its modulus 
a square of Lorentzian. Having the distribution function 
f(E = |E|) of local electric fields E in the domain pinning 
sites, we can derive the desired distribution function of 
switching times using the relation ((T|) . The same expres- 
sion ([T]) with respect to the fact that f(E) depends only 
on modulus E, which is a positive quantity, permits to 
determine uniquely the bounds of switching times spec- 
trum 



A/-Bcx 



(3) 



The expression 
Mcrtz's 



law 



for upper bound is consistent with 
131 ] where the characteristic time of FE 
polarization switching has the form exp(A/-E e xt) with 
A being activation field. In our approach parameter 
A = lIEo/hsT receives theoretical explanation. The ex- 
pression © also justifies theoretically Eq. (11) of Ref. 
[7j, its derivation will be given below. 

In the spirit of the statistical method [lH we write the 
field distribution function in the form 



/(E)=<5(E-E(n,..,r N )), 
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where bar means the averaging over spatial disorder and 
E(ri, tn) is the electric field created by N randomly 
distributed defects. At low pinning centers concentra- 
tion, where their spatial correlations (which means, for 
example, the absence of close pairs of defects) are small, 
the contributions of various defects to E are simply ad- 
ditive 
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E(r l! ...,r N ) = ^E(r i ). 



The dipole nature of defects implies that 



E(d,r) = -[3(l-n)n-l], 



(5) 
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where 1 = d/d and n = r/r arc unit vectors directed 
along d (vector of dipole moment) and r respectively. 
Using integral representation of 3D delta- function ((4]) and 
factorizing it by Eq. ([5]) with respect to (j6|), we obtain 
following expression for distribution function 
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Here n = N/V is pinning centers (dipoles) concentra- 
tion, d is its dipole moment. The distribution function 
([7| depends on vector E modulus due to spherical sym- 
metry of scalar products. The evaluation of integral ([7]) 



gives following explicit expression for field distribution 
function 



f(E) = 
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The expression © for f(E) along with Eq. ([T]) permits to 
obtain the explicit form of the switching times distribu- 
tion function g (In y) , y = r j t m . For that we demand that 
function g(hiy) should be normalized to unity, namely 

in ^ E\^)f{Em^P-d^ = 
Je(tP)=o dtp 

#)#=!, t[) = hiy = \ (9) 
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With respect to the fact that dE/dip < 0, Eq. © implies 
that ^ min = ip(E oo) = and Vmax = ip(E = 0) = 
X/E ext . This gives the above expression 

We emphasize here that for ordinary Lorentzian (we 
can obtain it supposing that electric field is a scalar quan- 
tity) the natural boundaries of switching times spectrum 
i"min and r max do not appear - both E (now it is a scalar) 
and lnr can be varied from — oo to oo. We also note 
that the case lnr = — oo is unphysical as it means zero 
switching time, i.e. system switches its polarization im- 
mediately. This is impossible as dipoles have finite mass 
so that some time should pass for polarization to reverse 
its direction. This means that consistent account for vec- 
tor nature of electric field gives physically correct switch- 
ing times spectrum in a ferroelectric. 

The explicit form of g(tp) can easily be extracted from 
Eq. © to read 
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The expression (flT)]) gives the desired distribution of 
switching times as a function of pinning centers (dipoles) 
concentration, temperature and external electric field. 
Note that distribution function (|10p depends also on the 
microscopic parameters of a material like pinning center 
dipole moment d as well as on parameters Eq and U of 
domain walls subsystem. 

The plots of g{ip) (in dimensionless variables A = 8/Eq 
= gjff, V = k B T/U, £ = E cxt /E ) are reported on 
Figff] at different fields, temperatures and A's. It is 
seen from FigfTJi (where r max and r m j n are also shown, 
which are constants in this case) that at small widths 
A the distribution function is peaked near r max , becom- 
ing delta function as A — > 0. This means that at small 
pinning centers concentration the systems switched non- 
randomly with the same time T max , which obeys Mcrtz's 
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law ([3]) [l3|. In the opposite limiting case the degree 
of disorder is maximal so that the distribution function 
is spread evenly between T m j n and r max . FigJTL shows 
the field dependence of the distribution function at fixed 
temperature and 5. According to the physics of our sys- 
tem, at high fields the distribution function is localized 
near = T mnl (note that at high fields r max — > T m - m , 
see Eq. ((21)), while at lower fields the distribution func- 
tion smears between = r m - m and r max which grows 
infinitely at field lowering. Thus, the electric field also 
effectively inhibits disorder in the system - at high fields 
the switching occurs at the same time Too ■ Note that this 
conclusion coincides with the results of Rcf. Q. Quali- 
tatively similar behavior is reported on FigJTJ;, where the 
temperature dependence of g(ip) (fT0| (at fixed field and 
S) is depicted. At hight temperatures the distribution 
function is localized near r m j n , while at low temperatures 
it becomes delocalized. This is because in the switch- 
ing process there is a competition between pinning and 
thermal energies of a dipole. At sufficiently high tem- 
peratures the thermal energy ksT exceeds the majority 
of energies of pinning centers so that the width of g(tp) 
in this case can be regarded as a measure of number of 
"turned on" pinning centers. Take it other way, at low 
temperatures, the random pinning in the system is strong 
and hence the distribution function is wide, while at high 
temperatures almost all pinning centers are "turned off" 
(their energies are smaller then the energies of dipoles 
thermal reorientations) and distribution function is nar- 
row. Thus, at sufficiently high fields and/or temperatures 
the system becomes effectively ordered and switching oc- 
curs coherently in time at r = = r m ; n . 

The expression (JTUJ) can be readily substituted into © 
to calculate switching polarization fraction p(t). Lat- 
ter calculation, however, can be done more conveniently 



in terms of field variable E. Namely, < e 



-t/r 
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4tt E 2 f(E)e- t/T{E) dE, where t(E) is given by Eq. 
((T]) . Proceeding along these lines we obtain following ex- 
pression for p(t) 
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where t\ = i/Too. Note that similar to Refs.0, EH we 
can pass to logarithmic variables and change exponen- 
tial function in the integrand pip to unit step function 
centered at zq = (1/(77 In ti) ~~ Q/A. This gives an ana- 
lytical approximation to the integral (|llj) (l5[. However, 
this analytical approximation is not always accurate. 

The results of numerical calculation of p(t) (JXTJ arc 
reported on Fig(5J On the same figure, as an illustrative 
example, we plot the experimental points from Ref. 0], 
where the switching kinetics in Pt/Pb(Zr,Ti)03/Pt thin 
films had been measured. It is seen, that distribution 
function (|10j). similarly to simple Lorentzian, provides a 
good coincidence between theory and experiment. As 
function ((TU)) decays at infinities like 1/ In 2 (t/t^) (this 
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FIG. 1: Distribution function of switching times as a function 
of: dimensionless width (product of dipoles concentration and 
its dipole moment) (a), electric field (b) and temperature (c). 
Note the common logarithms on the abscissa axis. 



decay character is similar to Lorentzian), the switching 
curve (|lip covers many decades in log(f) (here we con- 
sider the common logarithm of time). 

Our analysis shows, that expression (|lip for p(t) can 
(at different values of parameters 77 and £) describe 
also the experiments from Refs. [ijj Eal • This means that 
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FIG. 2: Comparison between our theoretical switching curve 
with experimental points from Ref. 0). 



17 ), augmented by relaxation term. Namely, the equa- 
tion of motion for the polarization vector P should have 
the form /xP + rjP = — 5F/8P, where F is correspond- 
ing (nonlinear) Ginzburg-Landau free energy including 
polarization gradients and, if necessary, the contribution 
from electrodes. This equation of motion should be aug- 
mented by the Maxwell equations for polarization outside 
the sample and proper boundary conditions accounting, 
in particular, for depolarization field. This problem is 
surely doable numerically, but for thin films it is doable 
analytically as domain structure there is always "sinu- 
soidal" [1S|] , i.e. it is described by the linear equation 
od state /(P) = — 8F/SP = 0. However, the proper ac- 
count for disordered ensemble of pinning centers seems 
to be problematic in such continuous theory. That is 
why the approaches based on the " improvement" of K AI 
model by the introduction of nucleus switching time dis- 
tribution seem to be reasonable for not perfectly ordered 
ferroelectric films. 



the above approach is at least not worse then previ- 
ous descriptions based on the phenomenological intro- 
duction of Lorcntzian distribution. For our fits, we take 
= 10~ 12 s, j] = 0.2 (the measurements in Ref.Q had 
been performed, most probably, at room temperature) 
and £ = 0.07, A = 0.07 for voltage v = 0.6 V, £ = 0.1, 
A = 0.2 for v = 1.35 V and £ = 0.3, A = 0.85 for v = 3 
V. The fact that we need larger A for fitting at higher 
fields is consistent with our theory as at increase of A the 
maximum of distribution function shifts towards Too, see 
FigOJ,. The higher fields, by themselves, make the sys- 
tem to switch coherently at smallest possible time t 00i 
which is seen on the Fig. [2j where the p(t) curve for 
?) = 3V has been switched already at the times less then 
10 -8 sec. The electric field may increase A via the dipole 
moments of pinning centers or their interaction with do- 
main walls. Theoretical explanation of this fact is out of 
frames of present paper. Actually, to have reliable base 
for creation of microscopic theory of domain nuclcation 
and switching in thin ferroelectric films, the detailed ex- 
perimental investigations of microscopic characteristics 
of domain walls motion are highly desirable. 

To conclude, we have shown that our model based on 
"squared Lorentzian" distribution function gives a good 
description of the switching kinetics in ferroelectrics thin 
films. Generally speaking, our treatment has to do with 
disordered ferroelectric films, where long-range random 
interaction is responsible for distribution of switching 
times. Even in ordered (epitaxial) films there are ran- 
domly situated technologically unavoidable defects (like 
stacking faults in PZT films, 0]), which could be the 
source of above distribution function of random fields 
and hence of switching times. We note here, that in 
perfectly ordered films the theory of ferroelectric switch- 
ing may be constructed on the base of continuous Chen- 
sky and Tarasenko approach [Tt} (see Eq. (38) of Ref. 
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